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Abstract 

We consider a directed percolation process on an M.xM rectangu- 
lar lattice whose vertical edges are directed upward with an occupation 
probability y and horizontal edges directed toward the right with oc- 
cupation probabilities x and 1 in alternate rows. We deduce a closed- 
form expression for the percolation probability P{x, y), the probability 
that one or more directed paths connect the lower-left and upper-right 
corner sites of the lattice. It is shown that P{x, y) is critical in the 
aspect ratio a = M-jM at a value ac{x,y) = [1 — y"^ — x(l — y)^]/2y^ 
where P{x, y) is discontinuous, and the critical exponent of the corre- 
lation length for a < acix,y) is ly = 2. 
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An outstanding unsolved problem in stochastic processes is the consider- 
ation of directed percolation [1, 2]. Directed percolation is a Markovian bond 
percolation process in which bonds are directed such that only clusters with a 
"flow" are relevant. Very few exact results of directed percolation are known. 
In 1981 Domany and Kinzel [3] solved one version of a directed percolation 
where the occupation probability is fixed at unity in one spatial direction of a 
rectangular lattice. The problem was subsequently reformulated and solved 
as a random walk by one of us and Stanley [4] . However, the Domany-Kinzel 
model is essentially of a one-dimensional nature due to the restricted freedom 
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in one spatial direction. To uncover the genuine nature of a two-dimensional 
directed percolation it is necessary to relax this uni-directional restriction. 

As a first step toward this goal we consider in this paper a directed perco- 
lation in which the unity percolation probability occurs in every other row of 
a rectangular lattice. We deduce a closed-form expression for the percolation 
probability and analyze its critical properties for large lattices. 

We first describe our model. Consider a 2-dimensional rectangular net of 
(M +1) X (2iV + l) sites with an aspect ratio 

a = M/2N. (1) 

Number the sites by (m, n) with m = 0, 1, • • • M, n = 0, 1, • • • 2A'" as shown 
in Fig. 1. Consider a bond percolation process on the lattice with vertical 
edges occupied with a probability Py — y and horizontal edges in the n-th 
row occupied with a probability 

Px — 1, n — odd 

— X, n — even. (2) 

Direct edges in the upward direction and toward the right. Occupied 
edges form directed paths if traced along the arrows. In ensuing discussions 
we shall refer to percolation configurations as bond configurations. A bond 
configuration is percolating if it contains one or more directed paths connect- 
ing the two opposite corner sites (0,0) and {M,2N). A typical percolating 
configuration is shown in Fig. 1. 

In a bond configuration there are {resp. MN — Ux) occupied (resp. 
empty) horizontal edges, and Uy {resp. 2(M -|- 1)A^ — Uy) occupied {resp. 
empty) vertical edges. Then the percolation probability, the probability that 
a bond configuration is percolating, is 

PMM^,y)= E a;-(l-a;)^^--y"^(l-y)2(^+^)^-"^ (3) 

perc conf 

where the summation is restricted to percolating bond configurations. It is 
clear that < PM,2N{x,y) < 1 since the summation (3) is identically 1 if 
unrestricted. It is also clear that Poo,2N{x,y) — 1 and Po^2N{x,y) — 0. Our 
interest is to investigate how does P change from 1 to as a varies, and 
whether the change is a sharp transition. 
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We state the main result as a Proposition: 
Proposition: 

For any x e [0, 1] and y G (0, 1), there exists a critical aspect ratio 

ac{x,y)^[l-y^-x{l-yf]/2y^ (4) 

such that 



lim P2aN,2N{x,y) = < 

Af— >oo 



' 1 if a> ac{x,y) 
if a < adx, y) (5) 
\ if ac{x,y). 



Moreover, for a < ac{x,y), we have the asymptotic behavior 

P{2aN,2N) e-^^/^ (6) 

where 

(7) 



Remarks: 

1. Equation (6) defines { as the correlation length and Eq. (7) gives the 
correlation length critical exponent 1^ — 2. 

2. For X — 1 our model reduces to the Domany-Kinzel model [3, 4] on an 
(M + 1) X (2A'" + 1) lattice and (4) leads to ac— {l — y)/y in agreement with 
previous result. For a; = our model is again a Domany-Kinzel model but 
on an (M + 1) x (A^ + 1) lattice with a vertical edge occupation probability 
y"^. Our result gives the critical aspect ratio 2ac = (1 — again in 
agreement with [3, 4]. 

Proof of the Proposition: 

The main body of this paper is the proof of the Proposition. 

There are 2N rows of vertical edges in the lattice. Number these rows 
from 1 to 2N starting from the bottom. An occupied vertical edge in a bond 
configuration is wet if it lies on a percolating path connecting (0, 0) and 
(M, 2A'") , and is primary wet if it is the first wet edge (in a row of vertical 
edges) counting from the left. In the bottom row of vertical edges in Fig. 1, 
for example, there are two wet edges and the primary wet edge is the one 
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connecting sites (1,0) and (1, 1). In a percolating configuration there is one 

primary wet edge in every row and these edges carry an overall occupation 
probability y'^^ . Since a bond configuration is percolating whenever a vertical 
edge in the 2A^-th row is primary wet, which can occur at any of the m-th 
horizontal positions m = 0, 1, • • • , M, we have 

M 

Pm,2n{x, y) = y^^ ^m,2iv- (8) 

m=0 

Here y'^^Wm,2n is the probability that the primary wet edge in the (2n)-th 
row occurs at the horizontal position m. 

We first establish a Lemma: 

Lemma: 

_ 1 /■ dt 

where the contour of integration is around the unit circle and 

a=l-y^ + x{l-yf, b = x{l-yf. (10) 

Proof of the Lemma: 

It is not difficult to see that the function Wm^2n{x, y) satisfies the recursion 
relation 

m 

Wm,2n = Ywk,2Wm~k,2n-2 (H) 
fc=0 

and the initial condition 

Wm,0^ ^Kr{m,0). (12) 

Define generating functions 

oo 

Wi{t) = E^--2r (13) 

m=0 
oo oo 

W2it,s) = EE^-,2ni"^«" (14) 



m=0 n=Q 



Substituting (11) into (14) and changing the order of summation by using 
J2m=o ^n=o — ^'n=o ^m=m obtain after some rearrangement and the use 
of (12), 

W2{t,s) = l + sWi{t)W2{t,s) 
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which yields 

W2(t,s) 

We can now invert (14) to obtain 



(15) 



'Wm,2n 





f dt 






(27ri)2 j 









1 

27n 



dt 



(16) 



where the contour of integration is around the unit circle. 

To compute Wi{t) we need to evaluate w„i2{x,y) for an (m + 1) x 3 
lattice. There are now 2 rows of vertical edges. As aforementioned y'^Wm,2 is 
the probability that (0, 0) is connected to (m, 2) with the primary wet vertical 
edge in the top row occurring at m. However the primary wet vertical edge 
in the bottom row can be at any j in < j < m. Denote the probability for 
this to occur by y^Xj{l — y)"^~^x"^~^. Then we have 



Wr 



.2 = E^i(l 

i=0 



(17) 



where the factor (1 — y)"^ ^ ensures that the primary wet edge in the 
top row is at m as shown in Fig 2(a). Particularly, we have wo,2 = Aq = 1. 

The factor Xj in (17) satisfies a recursive relation which can be written 



as 



m. 



= (1 - y)^j-i + 1/(1 - - y) j = 1, 2, 

The two terms on the right-hand side of (18) arise from the two possibilities 
that the vertical edge connecting (j — 1,0) and (j — 1,1) is either empty 
(with probability 1 — y) or occupied (with probability y) as shown in the two 
panels in Fig. 2(b). In the latter case the factor {1 — x){l — y) ensures that 
the site (j — 1, 1) is not on a percolating path. 

To solve the coupled recursion relations (17) and (18), define the gener- 
ating function 



A{t) 



j=Q 
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Multiplying (17) and (18) by and P ^, respectively, and summing over m 
and j — 1 from to oo, we obtain after some manipulation 



1 



t 

This gives 



l-xil-y)t' 

m-l] = {l-y)A{t) + y{l-x){l-y)W,{t). (19) 



- 13^ (20) 

after eliminating A(t) where a,b are given in (10). The substitution of (20) 
into (16) establishes the Lemma. 

We now continue the proof of the Proposition. 

Substitute (9) into (8) and carry out the summation in m. This leads to 
y^'^ f dt / 1 \ 

where the contour C+ is a circle enclosing the unit circle. Let ti and t2 
be the two roots of 1 — at + bt^ = 0, both of which are real. We have 
tit2 — i/b, ti + t2 — a/b, and hence 



{ti-l){t2-l) = tit2-{tl + t2) + l 

la y 
6~ 6^"^" fl-x^2 >0' 



so both ti and t2 lie outside the unit circle. We can therefore choose the 
radius of C+ to be greater than 1 but smaller than both ti and t2 so that 
C+ encloses only the simple pole t = 1 in (21). It follows that the first term 
on the right-hand side of (21) picks up only the residue at t — 1 which is 



y2N 



(l-a + b)^ 

and we obtain 

Pm,2n{x, y) = 1 — Im,n 
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where 



,2N 



dt 



^^'^ 2ni Jc+ {t - l)tA^+i(l -at + W^Y' ^^^^ 
Note that since |t| > 1 along C+ (22) leads to the expected result Poo,2N = 1- 
To further evaluate Im,n we introduce 2; = 1/t to write 



y 



2N 



(23) 



27ri Jc- (z - l)(z2 -az + b)^ 

where the contour C— is now within the unit circle. 

For M, N large and fixed aspect ratio a = M/2N, we can rewrite (23) as 



IM,N 



2iTi JC- z — 1 



dz r liV 
fa{z) 



(24) 



where 



fa{z) 



y^z^+'^ 
z'^ — az + b 



The integral Im,n can be evaluated using the method of steepest descent [5, 6] 
by deforming the contour to pass a point z = Zo where /^(-z) is stationary. 
To the leading order this gives Im,n ^ [/a(^o)]^- Moreover, since Im,n < 1, 
we must have fa{zo) < 1 with the equal sign holding at fa{zo) ~ 1. Thus a 
transition occurs at Zq — 1. 



Now 



OLZ^ — (1 + OL)a 2; + (2 + OL)b 



(^2 - + 6)2 
and the stationary point Zq is determined by 

az^ — (1 + OL)a 2^0 + (2 + OL)b — 0. 



The critical condition Zq = 1 now gives 

a -2b 



a — 



1-a + b 



(25) 



where ac{x, y) is given in (4). It is readily verified that we have {da / dz) z=i < 
along (25). Thus, for a > ac{x,y), the stationary point Zq lies within 
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the unit circle so we can deform C— continuously to pass Zq, and obtain 
haN,N = [fa{zo)]^ ~ 0. This givBS P2aN,N{x,y) ~ 1 which establishes the 
first line of (5). 

On the other hand, for a < ac{x,y), zq occurs outside the unit circle and 
when the contour C— is deformed to pass zq it must cross the simple pole at 
z — 1 and picks up the residue at the pole, which is equal to 1. This gives 
-^2aJV,Ar ~ 1 - [faizo)]^ and P2aN,N{x,y) ~ [/a(^o)]^ ~ for large N. This 
estabhshes the second hne of (5). 

For a — ac{x,y), zq is on the unit circle so the crossing of the contour 
at z = 1 picks up only half of the residue, namely, 1/2. This estabhshes the 
third line of (5). 

Finally, for a < ac{x,y), the method of steepest decent [5, 6] dictates 
that we have 

where expressions of Ci{x,y) and C2{x,y), which do not affect our conclu- 
sions, can be explicitly evaluated. This establishes the asymptotic behavior 
(6) with^ = 2/[C2(x,y)(a-ac)']. 

We have completed the proof of the Proposition. 

In summary, we have obtained a closed-form expression for the percola- 
tion probability Pm,2n{x, y) for the directed percolation process in which the 
occupation probability is y in the vertical direction and alternately x and 
1 in the horizontal direction. For M, N large, the percolation probability 
exhibits a critical behavior at a = ac- The correlation length ^ for a < ac 
is found to diverge with the critical exponent v = 2. While these properties 
are similar to those found in the Domany-Kinzel model [3, 4], our analysis 
permits the relaxation of the restriction of unit occupation probability in 
one spatial direction. It is hoped that the analysis serves as the first step 
of further relaxation in percolation probabilities, eventually leading to an 
understanding of genuine 2-dimensional directed percolation processes. 

This work was initiated while FYW was at the National Center of The- 
oretical Sciences (NCTS) in Taipei. The support of the NCTS is gratefully 
acknowledged. Work of LCC has been supported in part by the National 
Science Council, Taiwan. 



8 



References 



[1] R. Durrent, Oriented percolation in two dimensions, Annals of Probabil- 
ity 12, 999-1040 (1984). 

[2] B. D. Hughes, Random walks and random environments. Vol. 2 (Oxford 
University Press, New York 1996) p. 59. 

[3] E. Domany and W. Kinzel, Directed percolation in 2 dimensions: Nu- 
merical analysis and an exact solution, Phys. Rev. Lett. 47, 5-8 (1981). 

[4] F. Y. Wu and H. E. Stanley, Domany-Kinzel model of directed perco- 
lation: Formulation as a random-walk problem and some exact results, 
Phys. Rev. Lett. 48, 775-778 (1982). 

[5] P. M. Morse and H. Fcshbach, Methods of thoretical physics. Vol. 1 
(McGraw-Hill, New York 1953) p. 434. 

[6] P. Dennery and A. Krzywicki, Mathematics for physicists (Harper and 
Row, New York 1967), reprinted (Dover, New York 1996) p. 87. 

Figure Captions 

Fig. 1. A typical percolating configuration on a 6 x 5 lattice (M = 
5,A'" = 2). Open circles denote lattice sites. Oriented edges are occupied 
with weights shown. Empty edges carry weights 1 — x and 1 — y in horizontal 
and vertical directions respectively. 

Fig. 2. Construction of recursion relations, (a) Construction of (17). (b) 
Construction of (18). Occupied edges are shown as oriented edges; dotted 
edges can be either occupied or empty. To each row of vertical edges there 
is an additional factor y not shown. 
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Figure 1 




Figure 2 



